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Effect of Elastic Deformations on the Critical Behavior of Disordered Systems with 

Long-Range Interactions 
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Omsk State University, 55-a, pr. Mira, Omsk, Russia, 644077 \\ 
(Dated: February 2, 2008) 

A field-theoretic approach is applied to describe behavior of three-dimensional, weakly disor- 
dered, elastically isotropic, compressible systems with long-range interactions at various values of 
a long-range interaction parameter. Renormalization-group equations are analyzed in the two-loop 
approximation by using the PadBorel summation technique. The fixed points corresponding to 
critical and tricritical behavior of the systems are determined. Elastic deformations are shown to 
changes in critical and tricritical behavior of disordered compressible systems with long-range in- 
teractions. The critical exponents characterizing a system in the critical and tricritical regions are 
determined. 
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The effect of a long-range interaction described by a 
power law r DtT was studied analytically in the framework 
of the e-expansion [13] and numerically by Monte Carlo 
methods [46] for two- and one-dimensional systems. It 
was shown that effects due to long-range interaction are 
essential for the critical behavior of Ising systems when 
a < 2. The two- loop approximation applied in [7] di- 
rectly in the three-dimensional space corroborates the 
predictions of the £-expansion for homogeneous systems 
with long-range interactions. 

In structural phase transitions without piezoelectric ef- 
fect in the paraphase, elastic deformations play the role 
of a secondary order parameter whose fluctuations are 
not critical in most cases [8, 9]. Since the main contri- 
bution to striction effects in the critical region is due to 
the distance dependence of the exchange integral, only 
isotropic elastic systems are discussed below. 

It was shown in [10, 11] that coupling between order- 
parameter fluctuations and elastic deformations can be 
responsible both for change in critical behavior and for 
emergence of tricritical and tetracritical points. The in- 
troduction of frozen point impurities into a system both 
changes the critical behavior and eliminates multicriti- 
cal points [12]. The analysis presented in [13] showed 
that frozen structural defects introduced into spin sys- 
tems with long-range interactions increased the value of 
a long-range parameter corresponding to the transition 
to mean-field critical behavior. The effect of elastic defor- 
mations on homogeneous systems with long-range inter- 
actions also changes the critical behavior [14]. Therefore, 
it would be interesting to analyze the combined effect of 
structural defects and elastic deformations on systems 
with long-range interactions. 

In this work, the critical behavior of disordered com- 
pressible systems with long-range interactions in the 
three-dimensional space is described for various values 
of the long-range parameter a. 
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The Hamiltonian of the disordered Ising model includ- 
ing elastic deformations and long-range interactions can 
be written as 



+U a J d D qS q lS q 2S q 3S- q l^ q 2- q 3 

+a 3 J d D qy ql S q 2S- q i- q 2 (1) 
a (0) f 1 f 

+_ fr y ° / dDqS i s -v + 2 ai / dDq yi y ~i 

1 (0) r u 

1 di 9 iD i ^0 



where S q is an order parameter, it is a positive constant, 
To ~ \T — T c \/T c , T c is the phase-transition temperature, 
a is the long-range parameter, Ar q is a random impu- 
rity field (e.g., random temperature), a\ and a,2 are the 
elastic constants of a crystal, and a 3 is the quadratic 
striction constant. The coupling between impurities and 

3 

the nonfluctuating order parameter y(x) — ^ u aa {x), 

a=l 

where u a p is the strain tensor, is specified by the ran- 
dom field h q thermodynamically conjugate to u aa (x) In 
Eq. (1), integration is performed in the terms depend- 
ing on nonfluctuating variables that are not coupled to 
the S q , and the j/o terms (describing homogeneous de- 
formations) are separated. It was shown in [8] that this 
separation is necessary because inhomogeneous deforma- 
tions y q are responsible for exchange of acoustic phonons 
and for long-range effects that are absent under homo- 
geneous deformations. When the impurity concentration 
is low, random fields Ar g , h q and /iq can be treated as 



2 



Gaussian and specified by the function 

P[Ar,h, h ] = Acxpl--^- J Ar 2 d D q 

-kfW'-kf^ 9 
~kI ATqhqdDq ~ kl AT « h ° dDq] > 



result is 



(2) 



where A is a normalization factor and hi denotes posi- 
tive constants proportional to the concentration of frozen 
structural defects. 

Using the replica procedure to average over the random 
fields generated by frozen structural defects, we obtain an 
effective Hamiltonian of the system: 

. m 

A J a=l 

8 m f 

2 E / d D q(Sq 1 Sq 2 )(Sq 3 S b _ q i_ q2 - q3 ) 
a, 6—1 
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+.90 E / rfD <i^2^- 9 i- 9 2 (3) 

0=1"' 

JO) m r 1 r , , 

+J t^j dD ^ s % + \ x j d D qy q y- q + \^ 



where the positive constants S , go, 5o°\ ^ an d ^o can 
be expressed in terms of en and bi. The properties of 
the original system can be obtained in the limit m — > 0, 
where m is the number of replicas (images). 

Define an effective Hamiltonian depending only on the 
strongly fluctuating order parameter S q by the relation 

cxp{-H[S}} = B j exp{-H R [S,y}}l[dy q (4) 

If an experiment is carried out at constant volume, then 
Do is a constant and integration in Eq. (4) is performed 
over inhomogeneous deformations, whereas homogeneous 
deformations do not contribute to the effective Hamilto- 
nian. At constant pressure, the Hamiltonian is modified 
by adding the term Pfi, volume is represented in terms 
of the strain tensor components as 



fl = O [l + E U aa + E u aaU fj f) + 0(l 



5 )] (5) 



a=l 



and integration over homogeneous deformations is also 
performed in Eq. (4). It was noted in [15] that the 
quadratic terms in Eq. (5) can be important at high 
pressures and for crystals with large striction effects. The 
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The effective coupling constant v = u — z /2, which 
arises in the Hamiltonian owing to the striction effects 
determined by go, can have positive and negative val- 
ues. Therefore, this Hamiltonian describes both first and 
second-order phase transitions. When v = 0, the sys- 
tem exhibits tricritical behavior. Furthermore, the ef- 
fective interaction determined by the difference z w in 
Eq. (6) can also change the order of the phase transition. 
This representation of the effective Hamiltonian entails 
the existence of a higher order critical point where tricrit- 
ical curves intersect when vo = and zo — wq [16]. Note 
that Hamiltonian (6) is isomorphic to the Hamiltonian of 
the disordered Ising model with long-range interactions 
under the tricritical condition z = w - 

Using the standard renormalization-group procedure 
based on the Feynman diagram technique [17, 18] with 
the propagator G(k) = 1/(t + |fc| cr ), we derive the fol- 
lowing expressions for the functions P v ,Ps, Pz,Pw,lip and 
-ft specifying the differential renormalization-group equa- 
tion: 



-(2a -D)v 



1 - 36vJ + 245 J 
+1728 (2 Ji - J 2 - ^G)v 2 
-2304(2Ji - J 2 - ^G)vS 
+672 (2 Ji - J 2 - \g)5 2 
/3 S = -{2a- D)5 1 - 24vJ + 16<5J 
+576 (2 Ji - Jl - hj^v 2 - 

-1152(2 Ji - Jl - ^G)v5 + 352(2Ji - J 2 - 7^G)5 2 
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Z = -(2cr - D)z 1 - 24?; J - 2zJ + 86 Jq 

+576(2J! - J 2 - ^G)v 2 - 

-120(2Ji - J 2 - ?G)v<5 + 96(2Jx - J 2 - ^G)(5 2 
5 3 
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1 - 24v J + 8(5 J - 4z J + 2w J 
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In terms of the new effective interaction vertexes 



v 1 =v-J , v 2 = 6-J v 3 = z-J , v 4 = w-J . (8) 



the functions pi, j v and j t are expressed as 



/? 4 = -(2cr - £>)u 4 1 - 2Av 1 + 8v 2 - 4w 3 + 2v 4 
+576(2^ - 1 - \g)v\ - 120(2Ji - 1 - -G) Vl v 2 
+96(2Ji - 1 - lG)v 2 
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+288 (2 Ji - 1 - |Gjw 2 - 192(2 Ji - 1 - ^G) Vl v 2 

+32(2Ji-l-iG)« 2 , 
7v = (2(7 - D)64G(3i; 2 - 3wi« 2 + v 2 ). 



This redefinition is meaningful for a < D/2. In this 
case, Jo, Ji and G are divergent functions. Introducing 
the cutoff parameter A, we obtain finite expressions for 
the ratios J1/J5 2 and Gj Jq as A — > 00. 

The integrals are performed numerically. For a < D/2, 
a sequence of J1/J5 2 G/J5 2 corresponding to various val- 
ues of A is calculated and extrapolated to infinity. 

It is well known that perturbation-theory series are 
asymptotic and expressions (9) cannot be applied di- 
rectly since the interaction vertexes for fluctuations of 
order parameters in the fluctuation region are too large. 
For this reason, the required physical information was 
extracted from these expressions by applying the Pade- 
Borel method extended to the four-parameter case. The 
appropriate direct and inverse Borel transforms have the 
form 



f(v,6,z,w)= 



tl,...,U 



= fe t F(vit,v 2 t,v 3 t,V4t)dt, 


F(v,6,z,w) = £ + 1 .'.'."+t 4 )! t/ i U 2 ^3 "4 • 

To obtain an analytic continuation of the Borel trans- 
form of a function, a series in an auxiliary variable 9 is 
introduced: 



Pi = -(2cr - D) 
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(9) 
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k=0 



k\ 



(10) 

The Pade approximant [L/M] is applied to this series 
at 6 = 1. This procedure was proposed and tested in 
[19-22] for describing the critical behavior of a number 
of systems characterized by several interaction vertexes 
for order-parameter fluctuations. It was found in [19- 
22] that Pade approximants in the variable 9 preserve 
the symmetry of the system. This property is important 
for analysis of multivertex models. The approximants 
[2/1] are used to calculate the (3 functions in the two- 
loop approximation. 

Critical behavior is completely determined by the sta- 
ble fixed points of the renormalization group transfor- 
mation. These points can be found from the condition 



p i (v* 1 ,v* 2 ,v* 3 ,vt) = (i = 1,2,3,4). 



(11) 
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The requirement of stability of a fixed point reduces to 
the condition that the eigenvalues bi of the matrix 

dPi(v{,vZ,vt,v\) 
= • (12) 

are positive. 

The exponent v characterizing the growth of the cor- 
relation radius near the critical point (R c ~ \T — T c \~ u ) 
is determined by the relation 

a 

The Fisher exponent r\ describing the behavior of the 
correlation function near the critical point in the wave- 
vector space (G ~ k 2+v ) is determined by the scaling 
function y^: 77 = 2 — a + j v (vl, v%, v%, v%). The remaining 
critical exponents can be determined from the scaling 
relations. 

The table shows the stable fixed points of the renormal- 
ization group transformation and the fixed-point eigen- 
values of the stability matrix for a = 1.8 and 1.9. It was 
shown in [13] that stable fixed points exist in the physical 
region (v* > 0) for disordered systems only when a > 1.8. 
When a < 1.8, the stable points of any three-dimensional 
impurity system are characterized by negative values of 
the vertex vi. 



An analysis of the critical points and their stability 
shows that the fixed critical points of disordered systems 
with long-range interactions (nos. 1 and 4) are unsta- 
ble under elastic deformations. The critical behavior of 
disordered compressible systems with long-range interac- 
tions is described by their respective fixed points (nos. 2 
and 5) . The phase diagram of the substance can contain 
tricritical points specified by fixed points (nos. 3 and 6). 

Calculations of the critical exponents for the stable 
fixed points (nos. 2 and 5) yield 

cr = 1.9, v = 0.721, r} = 0.134 (13) 
cr = 1.8, v = 0.758, 77 = 0.251 

For the tricritical points (nos. 3 and 6), the critical ex- 
ponents are 

(7 = 1.9, v = 0.686, 77 = 0.134 (14) 
(7 = 1.8, v = 0.721, 77 = 0.251 

Thus, elastic deformations change both critical and tri- 
critical behavior of three-dimensional disordered Ising 
systems with long-range interactions. 

The work is supported by Russian Foundation for Ba- 
sic Research N 04-02-16002. 
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V*3 
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bi 


b 2 




bi 


a= 1.8 


1 


0.064189 


0.046878 








0.626* 


0.626* 


- 0.123 


- 0.123 
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0.064189 


0.046878 


0.066101 





0.626* 


0.626* 


0.124 


0.125 
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0.064189 


0.046878 


0.066101 


0.066101 


0.626* 


0.626* 


0.124 


- 0.124 


a = 1.9 
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0.066557 


0.040818 








0.559* 


0.559* 


- 0.118 


- 0.118 
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0.066557 


0.040818 


0.065716 





0.559* 


0.559* 


0.119 


0.119 
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0.066557 


0.040818 


0.065716 


0.065716 


0.559* 


0.559* 


0.119 


- 0.119 



